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ABSTRACT 

We have mapped the thermal emission line of SiO (v = 0; J = 2-1) associated 
with the quadrupolar molecular outflow driven by the very cold far-infrared 
source IRAS 16293—2422. The SiO emission is significantly enhanced in the 
northeastern red lobe and at the position ~50" east of the IRAS source. Strong 
SiO emission observed at ~50" east of the IRAS source presents evidence 
for a dynamical interaction between a part of the eastern blue lobe and the 
dense ambient gas condensation, however, such an interaction is unlikely to be 
responsible for producing the quadrupolar morphology. The SiO emission in the 
northeastern red lobe shows the spatial and velocity structure similar to those 
of the CO outflow, suggesting that the SiO emission comes from the molecular 
outflow in the northeastern red lobe itself. The observed velocity structure 
is reproduced by a simple spatio-kinematic model of bow shock with a shock 
velocity of 19-24 km s _1 inclined by 30-45° from the plane of the sky. This 
implies that the northeastern red lobe is independent of the eastern blue lobe 
and that the quadrupolar structure is due to two separate bipolar outflows. 
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The SiO emission observed in the western red lobe has a broad pedestal 
shape with low intensity. Unlike the SiO emission in the northeastern red lobe, 
the spatial extent of the SiO emission in the western red lobe is restricted to its 
central region. The spatial and velocity structures and the line profiles suggest 
that three different types of the SiO emission are observed in this outflow; the 
SiO emission arises from the interface between the outflowing gas and the dense 
ambient gas clump, the SiO emission coming from the outflow lobe itself, and 
the broad SiO emission with low intensity observed at the central region of the 
outflow lobe. 

Subject headings: ISM: individual (IRAS 16293-2422) - ISM: jets and outflows 
- ISM: molecules - stars: formation - shock waves 



1. INTRODUCTION 

Thermal emission from the SiO molecule is observed in shocked regions associated with 
molecular outflows from young stellar objects, in which the SiO abundance is enhanced by 
several orders of magnitude (e.g. Bachiller, Mart in- Pint ado, & Fuente 1991; Mikami et al. 
1992; Mart in- Pint ado, Bachiller, & Fuente 1992). Such an abundance enhancement of SiO 
is interpreted as the release of Si or SiO from dust grains by means of strong shocks (e.g. 
Mikami et al. 1992; Martin-Pintado et al. 1992). Since the SiO emission has been seldom 
detected in quiescent dark clouds because of its very low abundance (of the order of 1CT 12 ; 
Ziurys, Friberg, & Irvine 1989; Martin-Pintado et al. 1992), the SiO emission is generally 
believed to be a good tracer of shocked molecular gas. 

Strong SiO emission observed in low-mass star forming regions is mainly associated 
with outflows driven by very cold far-infrared sources referred to as the "Class 0" sources 
(e.g., Bachiller & Gomez- Gonzalez 1992; Bachiller 1996). On the basis of high-spatial 
resolution mappings of the SiO lines, it has been revealed that both spatial and velocity 
structures of the shocked gas are closely related to the history and an ejection mechanism 
of the outflows. The representative cases are the outflows driven by the Class sources in 
the LI 157 and L1448 dark clouds; the SiO maps present a picture that these outflows are 
driven by collimated and episodic jets ejected from the central star (Zhang et al. 1995; 
Gueth, Guilloteau, & Bachiller 1998; Dutrey, Guilloteau, & Bachiller 1997). Therefore, 
mapping of the SiO emission is expected to provide an opportunity to study the origin of 
the molecular outflows with complex morphology. 
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In this paper, we present the map of the SiO (v = 0;J = 2-1) emission in the 
quadrupolar outflow associated with IRAS 16293—2422 (hereafter referred to as 116293). 
116293 is known to be a proto-binary system with a projected separation of 840 AU (Mundy 
et al. 1992; Walker, Carlstrom, & Bieging 1993) embedded in the nearby (160 pc) p 
Ophiuchi molecular cloud. This outflow consists of two pairs of lobes; the brighter one 
extending along the east-west direction and the other one along the northeast-southwest 
direction (Walker et al. 1988; Mizuno et al. 1990). To explain the conspicuous quadrupolar 
morphology of this outflow system, Walker et al. (1993) proposed that the two pairs of lobes 
correspond to the two independent bipolar outflows driven by two independent sources. On 
the other hand, Mizuno et al. (1990) showed that the outflow is dynamically interacting 
with the dense ambient gas clump, and suggested that a single outflow lobe could be split 
into two lobes by the interaction. To examine which idea is plausible, and to probe the 
effect of the complex outflow on its surroundings, we investigate the spatial and kinematic 
structure of the shocked gas traced by the SiO emission. 



2. OBSERVATIONS 

Observations of the SiO (v=0; J— 2-1; 86.8 GHz) line were carried out with the 
Nobeyama 45 m telescope in 1992 April, 1993 April, and 1994 April. We used the SIS 
receiver with a single-sideband (SSB) system temperature of ~350 K. The half-power beam 
width and the main beam efficiency were 20" and 0.5, respectively. The backend was a 
bank of eight high-resolution acousto-optical spectrometers with a frequency resolution of 
37 kHz and a band width of 40 MHz, which correspond to a velocity resolution of 0.13 km 
s _1 and a velocity coverage of 140 km s _1 at 86.8 GHz. In addition to the SiO line, we also 
observed the H 13 CO + (J=l-0) line simultaneously with the same receiver. An area of 4.7' 
x 3.0' in Right Ascension and Declination was mapped with grid spacings of 14.1" and 
20". The map center was taken to be the position of the southeastern millimeter continuum 
source, a(1950) = 16 h 29 m 21.1 s and 5(1950) = -24° 22' 15.9". Offsets from this position 
are given as (Act, AS) in units of arcsecond. All observations were made in the position 
switching mode. Linear baselines were subtracted from the spectra. A typical rms noise 
level was ~0.16 K in T^. The telescope pointing accuracy was checked by observing the 
nearby SiO maser source toward V446 Oph at 43 GHz and was better than 5". 



3. RESULTS 
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3.1. Spatial Distribution and Kinematics of the Shocked Gas 

Figure 1 shows a map of the SiO emission integrated over Vlsr = 2 to 20 km s _1 
superposed on the CO J— 1-0 outflow map presented by Mizuno et al. (1990). The SiO 
J=2-l emission is detected in this velocity range. The SiO emission shows the U-shaped 
distribution open to the north with significant intensity enhancement in two clumps labeled 
El and E2 in Figure 1. The northeastern peak E2 (+90, +40) coincides with the peak of the 
northeastern red (hereafter referred to as NER) lobe in the CO map. On the other hand, 
the El peak (+50, 0) has no counterpart in the CO outflow map. The SiO emission is also 
detected at the position of 116293; however, its intensity is rather weak. At this position, 
detection of the higher transition lines of SiO J = 3-2, 5-4, 6-5, and 8-7 were reported 
by Blake et al. (1994) and Ceccarelli et al. (2000). The SiO emission observed toward the 
western red (WR) lobe shows two local intensity maxima. The position of one peak labeled 
Wl coincides with the peak of the high-velocity CO emission in the WR lobe, while the 
other peak W2 is located at the northern boundary of the WR lobe. No significant SiO 
emission is detected with Vlsr smaller than 2 km s -1 in the whole mapped region. Since 
most of the SiO emission is redshifted with respect to the cloud systemic velocity of ~ 4 km 
s _1 , the dense ambient gas interacting with the outflow is likely to exist behind the outflow 
lobes. 

Figure 2 shows a set of velocity channel maps of the SiO J = 2-1 emission from Vlsr 
= 2 to 20 km s _1 . It is shown that the El component is seen in the velocity range from 
Vlsr — 2 km s _1 to 10 km s _1 , while the E2 component contains higher velocity gas with 
maximum velocity reaching Vlsr — 20 km s™ 1 . Figures 2e-2i show that the peak position 
of the E2 component in the velocity range from Vlsr — 10 km s -1 to 20 km s -1 moves to 
the downstream as the velocity increases. 

The SiO line profiles observed at six representative positions, which are the IRAS 
position, El, upstream and downstream of E2, Wl, and W2, are plotted in Figure 3. The 
SiO line at the IRAS position (Figure 3a) peaks at Vlsr — 4.4 km s _1 with its main-beam 
brightness temperature of 0.56 K. The peak velocity of the SiO line is nearly the same as 
that of the H 13 CO + line. However, the line width of the SiO emission (AVfwhm = 3.5 km 
s _1 , where AVfwhm is the full width of half maximum) is much broader than that of the 
H 13 CO + emission (AV F whm = 11 km s -1 ), indicating that the SiO emission arises from the 
perturbed gas component, and not from the quiescent gas traced by the H 13 CO + line. On 
the basis of analyses using J = 5-4, 6-5, and 8-7 emission lines of SiO and its isotope, Blake 
et al. (1994) and Ceccarelli et al. (2000) have revealed that the SiO emission observed at 
the IRAS position arises from hot (T >60 K) and spatially compact (<3") region. A typical 
line profile observed in the El clump (Figure 3b) has a peak at the cloud systemic velocity 
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of Vlsr ~ 4 km s _1 with a gradual decrease in intensity toward the redshifted velocity and 
a rapid decline toward the blueshifted velocity. A line profile observed at the upstream of 
the E2 peak (+80, +40; Figure 3c) shows a double peak structure at Vlsr — 5 and 12 km 
s -1 , while that at the downstream of the E2 peak (+90, +70; Figure 3d) has its intensity 
maximum at Vlsr ~ 15 km s -1 . At the position of Wl (—20, 0; Figure 3e), the SiO line is 
pedestal shape without significant peak. The SiO emission at Wl is seen from Vlsr — 2 km 
s _1 to 13 km s _1 with typical brightness temperatures of 0.2 — 0.4 K. On the other hand, 
the SiO line observed near the W2 peak (—20, 60; Figure 3f) has a peak at the ambient 
cloud velocity with rather narrow line width of AVfwhm ~2 km s _1 . 



3.2. Comparison between the SiO and H 13 CO + distributions 

In Figure 4, we compare the spatial distributions of the low-velocity (Vlsr = 2-8 km 
s _1 ) and the high- velocity (V L sr = 8-20 km s" 1 ) SiO emission with the quiescent dense 
gas traced by the H 13 CO + emission (Vlsr = 2-6 km s _1 ). The H 13 CO + map shows that 
there are two dense gas condensations in the mapped area. One is associated with 116293 
and its peak is at (+10, —10). Toward this condensation, the velocity dispersion of the 
H 13 CO + emission AVfwhm is as large as ~1.4 km s _1 , which is the maximum value in the 
observed area. Another one is located at ~ 80" southeast of 116293. The H 13 CO + line width 
in this condensation is ~0.7-0.8 km s -1 , which is comparable to the typical line width of 
this molecule in the mapped region. This southeastern condensation is also observed with 
the NH 3 (1, 1) line by Mizuno et al. (1990). Hereafter we refer to this condensation as 
"16293E". 

Figure 4 shows that the spatial distribution of the SiO emission is anticorrelated with 
that of the H 13 CO + emission; especially, two prominent SiO clumps, El and E2, fall into 
the "valley" of the H 13 CO + emission. Since the spatial extent of the E2 clump is similar 
to that of the NER lobe, the intensity drop of the H 13 CO + toward the E2 clump suggests 
that the NER lobe is comprised of the ambient gas swept-up by the wind blown out from 
the protostar. It should be noted that the El clump that is dominated by the low-velocity 
emission is located between the two H 13 CO + condensations. Such a location of the El 
clump suggests that the SiO emission in this clump arises from the region in which the 
outflow driven by 116293 is interacting with 16293E. At the position of the El peak, the 
H 13 CO + line width has the second largest value of AVfwhm = 11 km s _1 . Such a line 
broadening supports the presence of dynamical interaction between the outflow and 16293E. 
The interaction between the outflow and 16293E is discussed in the next section. 
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4. DISCUSSION 

4.1. Dynamical Interaction between the Outflow and the Dense Gas Clump 

As mentioned in the previous section, the SiO emission in the El clump is likely to 
arise from the region where the outflow blown out from the protostar impinge on the 
dense ambient gas clump 16293E. Since the line profile observed in the El clump that 
shows a strong cut-off at the ambient velocity and a gradual redshifted wing (Figure 3b), 
it is suggested that the SiO emission in the El clump arises from the quiescent material 
accelerated by the shock. The dynamical interaction of the outflow with 16293E was also 
pointed out by Mizuno et al. (1990). They suggested that the blueshifted gas blown 
eastwards from its driving source is interacting with 16293E, and that such an interaction 
splits a single blueshifted lobe into the NER lobe and the eastern blue (EB) lobe, resulting 
in the quadrupolar morphology. However, the observed SiO emission is not blueshifted but 
redshifted; this clearly contradicts the above scenario. 

The redshifted SiO emission in the El clump suggests a possibility that the redshifted 
gas in the NER lobe is blown out to the eastward from its driving source and changes its 
direction to the northeast by means of a dynamical interaction with 16293E. In this case, 
the velocity field of the dense gas in 16293E is expected to be disturbed by the momentum 
provided by the outflow. On the basis of the CO observations, Mizuno et al. (1990) 
estimated the momentum of the outflowing gas in the NER lobe to be 0.5 M & km s^ 1 . 
Since the mass of 16293E is estimated to be ~1.6 M© from the H 13 CO + data by assuming 
the local thermodynamical equilibrium (LTE) condition with the excitation temperature 
T cx = 15 K and optically-thin H 13 CO + emission, the gas of 16293E is expected to show the 
velocity shift of ~0.3 km s" 1 toward the red from the rest of the cloud. However, there is 
no evidence of such a velocity shift in the H 13 CO + data. Therefore, we consider that the 
change of outflow direction by means of a dynamical interaction is unlikely. In addition, 
the velocity field observed in the E2 clump is explained by the bow shock impinging on the 
ambient medium (see §4.2), suggesting that the outflowing gas in the NER lobe is blown 
out to the northeast direction from its origin and does not change its direction on its way 
to the E2 position. 

Then, which outflow component is interacting with 16293E and is responsible for the 
strong SiO emission at the El position? A recent interferometric CO (J = 1-0) map has 
shown that the east-west pair of lobes, which corresponds to the EB and WR lobe pair, 
has an axis close to the plane of the sky and consists of fan-shaped lobes, each of which 
contains both blueshifted and redshifted components (Hirano et al. in preparation). This 
implies that the EB lobe contains a considerable amount of redshifted gas in the far side of 
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the lobe. Therefore, we consider that the redshifted gas in the EB lobe is interacting with 
16293E and is producing the redshifted SiO emission observed at El. 

4.2. Simple Kinematic Model of the NER Lobe 

As shown in Figure 1, the spatial distribution of the SiO clump E2 resembles that of 
the NER lobe in the single-dish CO maps. In particular, the map of the high-velocity SiO 
emission is quite similar to that of the high-velocity CO J = 2-1 emission presented by 
Ceccarelli et al. (1998). Channel maps presented in Figure 2 reveal that the high-velocity 
SiO emission in the E2 clump increases its velocity downstream. Such a velocity increase 
downstream is also observed in the CO (J = 2-1) emission in the NER lobe (Walker et 
al. 1988). These results indicate that the high- velocity SiO emission in the E2 clump has 
both spatially and kinematically similar characteristics to the CO outflow in the NER lobe, 
suggesting that this SiO emission comes from the NER lobe itself. Unlike the CO emission, 
the SiO emission in the E2 clump has another velocity component whose velocity is close 
to the cloud systemic velocity. The presence of two velocity components is explained if the 
low- and the high-velocity SiO emission components come from the near side and the far 
side of the NER lobe, respectively. Such two velocity components corresponding to the near 
sides and the far sides of the lobes were also seen in the position-velocity diagrams of the 
CO emission in the Mon R2 outflow (Meyers- Rice & Lada 1991). In the case of the NER 
lobe in the 116293 outflow, the low-velocity component has not been identified in the CO 
maps because it may be screened by overwhelming emission from the quiescent ambient gas. 

In Figure 5, we show the position- velocity diagram of the SiO emission along the major 
axis of the NER lobe (P. A. = 60°). The SiO emission in the NER lobe increases its velocity 
dispersion as an increasing distance from 116293, and the largest velocity dispersion is 
observed toward slightly downstream of the E2 peak. As pointed out by Bence, Richer, & 
Padman (1996), such a velocity feature is explained if the SiO emission arises from the bow 
shock generated by the high-velocity jet impinging onto the ambient medium. Recently, 
several authors have proposed the models that the molecular outflows are created through 
the propagation of the large bow shocks (e.g. Raga & Cabrit 1993; Masson & Chernin 
1993). In the following, we use a simple kinematical model of the bow shock proposed by 
Hartigan, Raymond, & Hartmann (1987) and compare the velocity features derived from 
the model with those observed. 

The geometry of the bow shock is schematically depicted in Figure 6. We assume 
that the bow shock has axial symmetry around the jet axis z with a shape described by 
a function z(r) = Ar 2 + _Br 4 , where r is the distance from the jet axis, and A and B 
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are numerical constants. In the frame of the exciting source, shock is propagating at a 
velocity of V s into the ambient pre-shock medium moving at Kmb- In the frame of the bow 
shock, material impacts the shock surface with velocity of V s — V am ^. Since the velocity 
component normal to the surface of the bow shock v n = (V s — l / am b)cos^ diminishes because 
of quick energy release via radiation, only the velocity component parallel to the surface 
f p = (V s — \4 m b)sin£ is conserved. Therefore, the post shock velocity is v 2 = (V s — \4mb)sin£, 
and its angle is 9 = (n/2) — £. The radial velocity V T in the frame of the exciting source is 
given by 

V T = (V s - Kmb)si<sin(^±0) - V s coscf), (1) 

where is the angle between the line of sight and the jet axis, and the + and — signs refer 
to the top and bottom halves of the curve in Figure 6, respectively. As demonstrated by 
Hartigan et al. (1987), the total velocity dispersion AV of the bow-shock feature and the 
median velocity Vidian are given by AV = V S - V amh and Vidian = -[(V 8 + V amh )/2] x cos0, 
respectively. In the case of the NER lobe, the total velocity dispersion is AV~16 km s^ 1 , 
and Kncdian measured in the frame of the exciting source (Vlsr ~4 km s _1 ) is ~8 km s _1 . 
To estimate the inclination angle 0, we compare the radial velocity distribution calculated 
from equation (1) with the observed position- velocity diagram. If we assume that the bow 
shock has a shape described as z(r) = 0.42r 2 + 0.136r 4 , which is designated as modified 
Raga shape in Hartigan et al. (1987), the observed feature is reproduced with 0~135° 
(i.e. the angle between the outflow axis and the plane of the sky % is ~45°), V s = 19.3 km 
s _1 , and V amb = 3.3 km s^ 1 (solid curve in Figure 5). In the case of the bow shock with a 
parabolic shape, z = r 2 , we obtained 0~12O° (i ~30°), = 24.0 km s _1 , and V am b = 8.0 
km s^ 1 (dash-dotted curve in Figure 5). Although the parameters depend on the assumed 
shape of the bow-shock surface, the observed velocity features can be explained if the axis 
of the flow is inclined by 30-45° from the plane of the sky and the shock velocity is ~19-24 
km s _1 . The dynamical time scale of the E2 shock is estimated to be (5-7) x 10 3 yr. It 
should be noted that the presence of the low-velocity component which is close to the cloud 
systemic velocity is a natural consequence in the context of the bow-shock outflow models. 

As shown in Figure 5, the velocity pattern observed in the NER lobe is reproduced by 
the simple bow-shock model, as in the case of the SiO emission observed toward the red 
lobe of the bipolar outflow in the L1448 dark cloud (Dutrey et al. 1997). This implies that 
the NER lobe is one of the bipolar outflow lobes and is not a splinter of the EB lobe. 

The bow-shock model cannot explain the high-velocity emission of Vlsr~18 km s _1 
observed farther downstream of the E2 peak (R.A. offset > 120"). Since the E2 peak is not 
located toward the end but toward the middle of the NER lobe, the bow shock at the E2 
peak is unlikely to be a terminal shock. Therefore, the NER lobe could have two-shock 
structure created by episodic ejection events as in the case of the blueshifted lobe of the 
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Ll 157 outflow (e.g. Zhang et al. 1995; Gueth et al. 1998). To assess this possibility, 
higher-resolution observations with interferometer are required. 

4.3. Three Types of SiO Emission Observed in the 116293 Outflow 

As discussed in the previous two subsections, the properties of the El and E2 clumps 
are different from each other. The SiO emission in the El clump arises from the interface 
between the outflow and the dense ambient gas clump, and that in the E2 clump is likely 
to come from the outflow lobe itself. 

The properties of the SiO emission observed at W2 appear to be similar to those of the 
emission from El, because the position of W2 corresponds to the periphery of the WR lobe 
at which the outflowing gas faces the ambient medium. The line profile of the SiO emission 
observed at W2 implies that most of the shocked gas has the velocity close to the ambient 
velocity. 

On the other hand, the SiO line profile observed at Wl, where the high- velocity CO 
emission is strongest in the WR lobe, shows a broad pedestal shape with its width measured 
at la level is 12 km s^ 1 . Unlike the NER lobe, the CO emission in the WR lobe shows no 
systematic displacement of emission peaks with increasing velocity, and has a peak at the 
position close to that of Wl at all velocities (Walker et al. 1988). The broad width of the 
SiO emission is likely to reflect the large velocity dispersion of the outflowing gas along the 
line of sight toward the Wl position. Unlike the E2 clump, the spatial extent of the SiO 
emitting region in the Wl is restricted to the small area of ~20" in size, which is much 
smaller than the size of the WR lobe observed in the CO. Such a broad SiO emission with 
low intensity is also observed in the central region of the outflows in S140 and IC1396N 
(Hirano et al. in preparation). 

To summarize, in the 116293 outflow, the SiO emission lines having three different 
properties are observed. The first one is the enhanced SiO emission at the interface between 
the outflowing gas and the dense ambient gas clumps (El and W2), which shows the line 
profile peaks at the velocity close to the cloud systemic velocity. The second one is the 
SiO emission showing similar spatial extent and kinematic structure as the CO outflow, 
and considered to arise from the outflow lobe itself (E2). The third one is the broad SiO 
emission with low intensity observed toward the central region of the outflow lobe (Wl). 
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4.4. SiO Abundances 

To estimate the SiO abundance in the shocked gas, we assumed optically thin emission, 
and calculated the column densities of the SiO molecules by assuming the LTE condition. 
The permanent dipole moment of the SiO molecule is employed to be 3.1 Debye. In the 
following, we assumed the excitation temperature T cx to be 15 K, which is almost the same 
as the rotational temperature estimated from the NH 3 (1,1) and (2, 2) transitions (Mizuno 
et al. 1990). It is possible that the excitation temperature of the SiO in the shocked gas 
is higher than that assumed here, because recent observations of high-excitation ammonia 
have revealed that the shocked gas in several outflows is heated to more than 100 K (e.g. 
Bachiller, Martin-Pintado, & Fuente 1993; Umemoto et al. 1999). If we adopt 100 K as 
an excitation temperature, the column density of the SiO becomes a factor of three higher 
than that calculated for T ex = 15 K. 

By using the total integrated intensities (from Vlsr = 2 to 20 km s" 1 ), the SiO column 
densities, iV(SiO), at the El and E2 positions are derived to be 4.3 xlO 13 cm -2 and 8.4 xlO 13 
cm~ 2 , respectively, whereas it is estimated to be ~ 2.0 x 10 13 cm -2 at both the Wl and 
W2 positions. At the IRAS position, we obtain A r (SiO) to be 1.0 x 10 13 cm -2 . When we 
adopt an H 2 column density averaged over the core of 0.2 pc in radius, 3.6 xlO 22 cm~ 2 , 
derived from the C 18 observations of Mizuno et al. (1990), the SiO abundance X(SiO) 
is estimated to be 1.2 xl0~ 9 toward the El peak, 5.2 xl0~ 10 toward Wl and W2, and 
2.8 xl0~ 10 at the IRAS position. It should be noted that the abundance values estimated 
here are the beam-averaged values. If the SiO emission observed at the IRAS position 
comes from the small area of <3" size, the X(SiO) becomes a factor of ~50 higher than 
that estimated here, which is consistent with the SiO abundance value derived by Ceccarelli 
et al. (2000) from their LVG analysis. 

To estimate the SiO abundance in the outflowing gas in the NER lobe, we compare the 
column density of the high- velocity SiO (from 8 to 20 km s _1 ) with that of H 2 derived from 
the CO data observed by Mizuno et al. (1990). At the E2 peak, the JV(SiO) and N(H 2 ) of 
the high- velocity gas are estimated to be 7.3 xlO 13 cm~ 2 and 8.8 xlO 20 cm -2 , respectively. 
We thus obtain X(SiO) in the high- velocity gas to be 8.3 xl0~ 8 at the E2 peak. This 
indicates that the SiO abundance in the NER lobe is comparable to that observed in the 
blueshifted lobe of the L1157 outflow (Mikami et al. 1992) and is enhanced by a factor of 
~10 4 as compared to the upper limits in the dark clouds (Ziurys et al. 1989). If we assume 
that all the Si atoms (Si/H 2 ~7xl0~ 5 ) is depleted on grains in the unperturbed component, 
at least ~0.1 % of Si in the interacting region should be released in the gas phase in order 
to explain the SiO abundance in the NER lobe. 
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5. Summary 

We have mapped the thermal emission line of SiO (J = 2-1) associated with the 
quadrupolar molecular outflow driven by 116293. Our main results are the following: 

1. The SiO emission is significantly enhanced in the northeastern red lobe (E2) and at 
the position ~50" east of the IRAS source (El). The SiO emission is also detected at the 
IRAS position and toward the western red lobe (Wl and W2), however, the line intensities 
observed at these positions are less than half of those observed at El and E2. 

2. Most of the SiO emission is redshifted with respect to the cloud systemic velocity 
of Vlsr = 4 km s _1 . This suggests that the dense ambient gas interacting with the outflow 
exists behind the outflow lobes. 

3. The El clump is dominated by the low-velocity emission (Vlsr = 2-10 km s _1 ) 
and is located between the dense gas condensation associated with 116293 and the other 
one 16293E, which is located at ~80" southeast of 116293. The location of the El clump 
and the SiO line profile observed there suggest that the SiO emission in the El clump arises 
from the region where the outflow blown out from 116293 is dynamically interacting with 
16293E. 

4. In the E2 clump, high velocity SiO emission with maximum velocity reaching Vlsr 
= 20 km s _1 is observed. The spatial and velocity structure of the SiO emission in the E2 
clump are similar to those of the NER lobe observed in the CO emission, suggesting that 
the SiO emission in the E2 clump comes from the NER lobe itself. The velocity pattern of 
the SiO in the NER lobe is reproduced with a simple spatio-kinematic model of bow shock 
with a shock velocity of 19-24 km s -1 inclined by 40-45° from the plane of the sky. This 
model gives a dynamical timescale of ^5000-7000 yr for the NER lobe, being consistent 
with the significant increase of the SiO abundance by a factor of ~10 4 relative to the values 
in dark clouds. 

5. Although the SiO emission at El presents evidence for a dynamical interaction 
between a part of the EB lobe and 16293E, such an interaction is unlikely to be responsible 
for producing the quadrupolar morphology. Since the velocity structure of the NER lobe, 
which is reproduced by the bow-shock model, is similar to that of the lobe in the bipolar 
outflow, the NER lobe is likely to be one of the lobes of the bipolar outflow rather than 

a splinter of the EB lobe. Therefore, we consider that contribution of two independent 
bipolar outflows is the plausible explanation for the origin of the quadrupolar structure of 
the 116293 outflow. 

6. In the 116293 outflow, the SiO emission lines with three different properties are 
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observed. They are 1) the enhanced SiO emission at the interface between the outflow and 
dense ambient gas (El and W2), 2) the SiO emission arises from the outflow lobe itself 
(E2), and 3) the broad SiO emission with low intensity observed toward the central region 
of the outflow lobe (Wl). 

We would like to thank the staff of NRO for the operation of our observations and 
their support in data reduction. N.H. acknowledges support from a Grant-in- Aid from the 
Ministry of Education, Science, Sports, and Culture (No. 09640315). 
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Fig. 1. — Integrated intensity map of the SiO J = 2-1 emission (Vlsr = 2-20 km s _1 ; gray 
scale) superposed on the CO J = 1-0 outflow map (Vlsr = —10-2 km s _1 ; blue contours 
and Vlsr — 6-20 km s -1 ; red contours) presented by Mizuno et al. (1990). The contours of 
the SiO J = 2-1 map are drawn every 1.0 K km s -1 (3 a) with the lowest contours of 1.0 
K km s" 1 (3 a). Two crosses indicate the positions of the 116293 binary protostars. Green 
line denotes the cut of the position- velocity diagram presented in Figure 5. 

Fig. 2. — Channel-velocity maps of the SiO J = 2-1 emission with the velocity interval of 2 
km s _1 . Contours are every 0.33 K km s _1 (3 a) with the lowest contours of 0.33 K km s _1 
(3 a). Solid and dashed grey lines denote the areas of the red and blue lobes of the CO J = 
1-0 outflow (Mizuno et al. 1990), respectively. The positions of two protostars are indicated 
with crosses. 

Fig. 3. — Line profiles of the SiO J = 2-1 emission observed toward (a) 116293, (b) El, (c) 
the upstream of the E2 peak, (d) the downstream of the E2 peak, e) Wl, and f) W2. The 
offsets are in arcseconds from the position of the southeastern continuum source: a(1950) = 
16 h 29 m 21. s l, 5(1950) = -24° 22' 15.9". 

Fig. 4. — Comparison of the high- velocity SiO (Vlsr = 8-20 km s -1 ; red contours), the 
low- velocity SiO (Vlsr — 2-8 km s -1 ; green contours), and the H 13 CO + J = 1-0 integrated 
intensity (Vlsr — 2-6 km s -1 ; gray scale) maps. The contours for the high- velocity SiO are 
drawn every 0.69 K km s _1 (3 a) from 0.69 K km s _1 (3 a), and those for the low- velocity 
SiO are every 0.60 K km s _1 (3 a) from 0.60 K km s -1 (3 a). The contours of the H 13 CO + 
map are drawn every 0.54 K km s" 1 (3 er) with the lowest contours of 0.54 K km s -1 (3 a). 
The crosses denote the positions of the protostars. 

Fig. 5. — Position-velocity diagram of the SiO J = 2-1 emission along the major axis 
of the NER lobe (PA = 60°; green line in Figure 1). The lowest contour level and the 
contour spacing are 0.3 K. The solid and dash-dotted curves represent the velocity patterns 
derived from the bow-shock outflow models with modified Raga shape and parabolic shape, 
respectively. 

Fig. 6. — Geometry of the bow shock impacting a downstream medium moving at V am b in 
the frame of the exciting source. The shape of the bow shock is given by z(r). <f) is the 
inclination angle between the axis of the bow shock and the line of sight (0 = n/2 + i). The 
shock propagation velocity is V s . 



